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ABSTRACT 

We present the results of a set of numerical simulations of long-duration gamma-ray burst jets associated with 
massive, compact stellar progenitors. The simulations extend to large radii and allow us to locate the region in 
which the peak frequency of the advected radiation is set before the radiation is released at the photosphere. 
Light curves and spectra are calculated for different viewing angles as well as different progenitor structures 
and jet properties. We find that the radiation released at the photosphere of matter-dominated jets is able to re- 
produce the observed Amati and energy-Lorentz factor correlations. Our simulations also predict a correlation 
between the burst energy and the radiative efficiency of the prompt phase, consistent with observations. 
Subject headings: ganmia-ray: bursts — hydrodynamics — methods: numerical — relativity 



1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) are characterized 
by bright, non-thermal radiation spectra with typical photon 
energies of a few hundred keV to several MeV (Band et al 
1993; Kaneko et al. 2006; Zhang et al. 2011). The stan- 
dard model to explain these observations invokes collision- 
less shocks between parts of the relativistic outflows moving 
at different speed (Rees & Meszaros 1994) to produce a tan- 
gled magnetic field and non-thermal electrons, eventually re- 
leasing synchrotron radiation (Daigne & Mochkovitch 1998; 
Piran 1999, Bosnjak et al. 2009). Alternatively, the mag- 
netic field may be advected from the central engine and the 
collisions working as a trigger for magnetic dissipation (the 
ICMART model, Zhang & Yan 201 1). 

When GRB observations are corrected for redshift and con- 
sidered in the burst's frame, a correlation between the burst's 
isotropic equivalent energy and its characteristic photon en- 
ergy is revealed (Amati et al. 2002). The correlation was orig- 
inally discovered with data from BeppoSAX, and dubbed the 
Amati correlation. Even though the role of selection effects 
has not been completely understood (Nakar & Piran 2005; 
Kocevski 2012), the correlation has been confirmed with data 
from all subsequent observatories (Amati 2006; Ghirlanda et 
al. 2008; Amati et al. 2009). To lend credibility to an intrin- 
sic origin, a correlation between burst peak luminosity and 
peak photon energy is also observed among bursts (Yonetoku 
et al. 2004; Ghirlanda et al. 2012a) and within individual 
events (Lu et al. 2012). The latter correlation is free of selec- 
tion effects and points to a physical origin for all the others. 
More recently, the analysis of a sample of GRBs with very 
early afterglow observations revealed a correlation between 
the isotropic equivalent energy of the prompt emission and the 
Lorentz factor of the outflow (Liang et al. 2010; Ghirlanda 
et al. 2012b). Taken together, these correlations present a 
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challenge to the standard synchrotron shock model (SSM). 
Outflows characterized by large Lorentz factors should pro- 
duce internal shocks at large radii, therefore generating mag- 
netic fields of much lower intensity, and eventually bursts 
characterized by smaller photon energies, predicting an anti- 
correlation in contrast to observations (Ghisellini et al. 2000). 

An alternative model for the prompt emission of GRBs is 
the photospheric model (Meszaros & Rees 2000; Pe'er et al. 
2005, 2006; Rees & Meszaros 2005; Giannios 2006; Lazzati 
et al. 2009). It does not specify how the photons are produced, 
concentrating instead on how the interaction of the radiation 
field with the leptonic component of the outflow modifies the 
spectrum in the optically thick phase, before it is released at 
the photosphere. The photospheric model has been proposed 
as a viable candidate to reproduce correlations such as the 
Amati correlation on the qualitative level (Thompson et al. 
2007; Lazzati et al. 2011; Fan et al. 2012). On the quanti- 
tative side, comparison of the burst properties from numeri- 
cal simulations to observations have shown that while photo- 
spheric emission from baryonic jets can reproduce the slope 
of the correlation, they cannot reproduce the normalization, 
simulated bursts being more energetic than observed ones for 
a given peak frequency (Lazzati et al. 2011; Nagakura et 
al. 2011). Deeper theoretical understanding of photospheric 
emission has revealed that radiation and matter are not in equi- 
librium all the way to the photosphere and therefore the peak 
of the spectrum is formed at moderate optical depths (r ~ 50) 
when the radiation is at higher temperature (Giannios 2012). 

In this paper we present results of an extensive set of sim- 
ulations aimed at exploring quantitatively the predictions of 
the photospheric model for baryonic jets. In contrast to our 
previous work, the peak frequency of the spectrum is calcu- 
lated at an optical depth larger than unity, following Giannios 
(2012). We also explore the effect of the jet injection proper- 
ties and of the progenitor structure on the ensuing light curves 
and spectra. This paper is organized as follows: in Sect. 2 
we present our numerical simulations, in Sect. 3 we detail the 
post processing methods to derive light curves and spectra, in 
Sect. 4 we present our results, and in Sect. 5 we summarize 
and discuss this work. 

2. NUMERICAL SIMULATIONS 

All the simulations presented in this paper were performed 
with the FLASH code (Fryxell et al. 2000), version 2.5. The 
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Fig. 1 . — Density profiles of the GRB progenitor stars used in this work. 

FLASH code is a modular block-structured adaptive mesh re- 
finement code, parallelized using the Message Passing Inter- 
face (MPI) library. It solves the Riemann problem using the 
Piecewise Parabolic Method (PPM). From the point of view 
of the research presented here, the main strength of this code 
is that it can perform special relativistic hydrodynamic com- 
putations on an adaptive mesh. This is crucial for the mod- 
eling of the interaction of the collimated jet with the stellar 
envelope. This approach allows for significantly enhanced 
resolution along the jet axis and the jet- star boundary with- 
out the need to resolve the rest of the star to the same extent, 
thereby making high-resolution simulations feasible. In order 
to increase the resolution near the center of the star and close 
to the axis of the jet we modified FLASH to be able to vary 
the maximum level of refinement allowed over the simulation 
grid (see also Lopez-Camara et al. 2013). Additional modifi- 
cations to the standard FLASH code are described in Morsony 
et al. (2007). 

We adopted a maximum resolution of 4 x 10^ cm at the 
highest level of refinement. At this resolution the transverse 
dimension of the injected jet is resolved into 44 elements. Our 
simulations do not include magnetic fields, due to the techni- 
cal challenge of performing MHD calculations with relativis- 
tic motions on an adaptive mesh. In addition, gravity from a 
central mass and self-gravity are neglected since the charac- 
teristic times of the jet- star interaction are much shorter than 
the dynamical time of the progenitor star's collapse. 

All our simulations adopted a realistic GRB stellar pro- 
genitor, and we explored three different progenitors all taken 
from Woosley & Heger (2006). Model 16TI is a 16 solar- 
mass Wolf-Rayet star with an initial metallicity 1% solar and 
angular momentum / = 3.3 x 10^^ erg-s. The mass of the 
star at pre-explosion is 13.95 solar masses and its radius is 
4.07 X 10^^ cm, corresponding to 0.6 solar radii. A rela- 
tivistic jet is injected as a boundary condition at a distance 
Rq = 10^ cm from the star's center. The jet has initial Lorentz 
factor Fo = 5 and an half-opening angle Oq = 10° (5° in one 
case). The jet is injected hot, with enough internal energy 
to allow for acceleration up to an asymptotic Lorentz fac- 
tor Too = 400 upon complete, non-dissipative acceleration 
(Foo = 100 upper limit was adopted in one case). The en- 
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Fig. 2. — False-color stratification map of the logarithm of the density of 
our fiducial 16TI simulation at f = 100 s. The thick black line shows the 
location of the photosphere for an observer located at a viewing angle Oo = 
0.5° towards the bottom of the panel. 
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16TI 
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100 
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350B 
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400 


10 


5.33 X 10^0 


100 
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160I 
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400 


10 


5.33 X 10^0 


100 
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16TI 
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400 


10 


10^0 


100 
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16TI 
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400 


5 


5.33 X 10^0 


100 
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16TI 


5 


400 


10 


5.33 X 10^0 


67 


8 


16TI 


5 


400 


10 


5.33 X 10^0 


30 



TABLE 1 

Details of the simulations presented in this work. 

gine luminosity is kept constant for a duration ^eng (typically 
^eng = 100 s, but shorter engines were explored as well), af- 
ter which the jet is turned off and the boundary condition is 
set to reflective. The bottom and polar boundaries are set to 
reflective for the whole duration of the simulation, while the 
two outer boundaries are set to absorbing. The simulation 
box is 2.5 X 10^^ cm in length (along the jet direction) and 
5 X 10^^ cm across, approximately a factor 10 larger than our 
previous simulations. 

Two additional progenitor stars were tested. Model 1601 is 
a star of 16 solar masses with an initial metallicity 10% solar. 
Its pre-explosion mass is 12.21 solar masses and its radius 
is 0.7 solar radii. Model 350B is a 35 solar-mass star with 
initial metallicity 10% solar, pre-explosion mass 21.24 solar 
masses, and radius 0.6 solar radii. The density profiles of the 
three progenitors are shown in Figure [T] while details of the 
simulations are reported in Table 

3. SPECTRA AND LIGHT CURVES 

The calculation of the light curves and spectra is per- 
formed analogously to Lazzati et al. (2011, see also Mizuta 
et al. 2011; Nagakura et al. 2011) with a few important 
changes. First, we changed the optical depth at which we 
compute the radiation temperature. In our previous publi- 
cations, we used the photospheric radius as the location at 
which the radiation temperature was calculated through the 
black body relation u = aT^, where u is the energy density 
and a = 7.56 X 10~^^ erg cm~^ is the radiation energy con- 
stant. We now use Eq. 6 from Giannios (2012) to derive the 
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opacity at which the radiation temperature is set: 



tt =46 



jl/6 .1/3 
^ 2.5 ^ '12.5 



(1) 



where f± is the number of leptons per proton (including 
electron-positron pairs), e is the fraction of the energy of the 
outflow carried by radiation, r] = L/Mc^ is the asymptotic 
Lorentz factor, and we used the notation = 10~^Q. Because 
of the decrease in the radius with respect to the standard as- 
sumption r = 1, the location at which the spectrum is formed 
is easier to determine from the simulation frames and no cor- 
rection needs to be applied (see Eq. 4 in Lazzati et al. 201 1). 
Our numerical simulations do not allow us to compute the 
number of electron-positron pairs and we assume that f± = l 
for the remainder of this paper. This assumption is justified by 
the fact that the comoving temperature is of only a few keV or 
less where the peak frequency is computed and, in absence of 
non-thermal radiation, pairs are not important. If a significant 
population of pairs were present, the location of the spectrum 
formation would be moved outwards and the peak frequency 
would be decreased, weakening the agreement with observa- 
tions (see figures and discussion below). 

The location of the region of the spectral peak formation 
was calculated analogously to Lazzati et al. (201 1) but updat- 
ing their condition to: 
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where P = /3(^iab,-^,z) is the local velocity of the outflow in 
units of the speed of light, T = r(^iab,x,z) is the local bulk 
Lorentz factor, and 6^ = Oy(t\^\)^x^z) is the angle between the 
velocity vector and the direction of the line of sight, x is the 
coordinate perpendicular to the line of sight, while z is the 
coordinate along the line of sight. All the values of /3, F, and 
are evaluated at the same delayed coordinate (^lab,-^,^) = 



(^obs ■ 



X, z) as the comoving density. 



Figure [2] shows an image of the density stratification from 
our fiducial simulation (simulation number 1 in Table [T]) at a 
time ^ = 100 s after the engine onset. A black line is used to 
show the location of the photosphere for an observer located 
at a viewing angle = 0.5° towards the bottom of the panel. 
Most of the photospheric emission comes from the bottom of 
the trough where the outflow velocity vector and the radiation 
propagation are aligned within an angle ~ 1/F. 

In addition, we refined the way in which we compute the 
peak frequency to make the comparison with observations 
more accurate. First, we consider the fact that what is mea- 
sured from the data is the peak of the uFiu) spectrum. Sec- 
ond, we consider the fact that the location of the peak fre- 
quency is affected by the addition of high-frequency power- 
law tails likely due to Comptonization (Pe'er et al. 2006; Gi- 
annios & Spruit 2007; Beloborodov 2010; Lazzati & Begel- 
man 2010) and low-frequency power-laws likely due to a syn- 
chrotron component (Vurm et al. 2011). The combination of 
these two effects shifts the peak to higher frequencies by ap- 
proximately a factor 2, depending on the slope of the high- 
and low- frequency power-law components (see Figure [3]). 

A sample of light curves, photospheric radii, and peaK fre- 
quencies are plotted vs. time in Figure]?] These light curves 




Fig. 3. — Modification of the location of the uFiu) peak of a Planck spec- 
trum due to the addition of a Comptonized high-frequency tail (blue) or low- 
frequency tail (red). Indices in the x-axis are spectral indices (F(iy) oc u'^; 
F{u) oc ly-^). 
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Fig. 4. — Sample results of our post processing for simulation #1 (see Ta- 
ble^. Each panel shows a different viewer. Different colors show the bolo- 
metric light curve (red), the radius of the photosphere (blue) and the peak 
frequency (green). 

are quite representative also of the other simulations. As a 
matter of fact, the effects of the viewing angle are much more 
pronounced than those of a different progenitor, engine dura- 
tion, or jet properties. As we move away from the jet sym- 
metry axis, all bursts grow weaker, start at a later time, are 
released at a larger photospheric radius, and are characterized 
by a smaller photon frequency. 

4. CORRELATIONS 

Figure |5] shows the observed Amati correlation (Amati et 
al. 2002; 2009) compared to the one found from our synthetic 
light curves. The change in the location of the radius at which 
the spectrum is formed has great influence on the synthetic 
light curves. The typical photon frequency is increased while 
the total energy is almost unchanged. As a consequence, syn- 
thetic bursts agree quantitatively with the Amati correlation 
in both the slope and the normalization. The agreement re- 
sults from the combination of two effects. First, the shearing 
and shocking effects of the stellar material on the outflowing 
plasma cause a polar stratification of the jet. Different ob- 
servers see different bursts that, when placed on the Amati 
diagram, form a stripe along the correlation. In addition, the 




Fig. 5. — Overlay of the simulation results on the observed spectral and 
energetic properties of gamma-ray bursts. The black symbols with error 
bars show the observed Amati correlation, i.e., the correlation between the 
isotropic-equivalent energies of the bursts and the peak photon energies of 
their uFiu) spectra (Amati et al. 2002; 2009). The colored lines and symbols 
show the results of our simulations. Different colors refer to different sets 
of progenitor and jet properties, as described in the upper left legend. The 
legend states explicitly the parameters of the simulation that differ from our 
fiducial progenitor/jet pair (see text and Table [TJ. Different symbols refer to 
different orientations of the observer with respect to the jet axis, as described 
in the lower right legend. 

pressure exerted by the stellar material on the jet walls col- 
limates the outflow in such a way as to produce large-scale 
jets that are insensitive to their initial conditions. We also find 
that the most important parameter to determine the location of 
a burst on the Amati diagram is the radius of its stellar progen- 
itor: the slightly bigger 1601 progenitor produces higher peak 
frequencies at a given energy than the more compact 16TI and 
350M progenitors. 

The agreement between the synthetic bursts and the ob- 
served correlation is quite remarkable since the photospheric 
model has only a few free parameters: the progenitor star 
structure, the outflow luminosity and its asymptotic Lorentz 
factor, and the opening angle of the injected jet. We tried 
different values for each of these parameters and we always 
found complete agreement between the simulated and ob- 
served bursts. We also tried bursts from engines with different 
durations. As long as the jet reaches the star surface still rela- 
tivistic, the computed bursts lie on the Amati correlation. We 
finally note that our simulations do not cover extreme bursts 
with very high luminosity and Lorentz factors but are rather 
targeted to reproduce "average" events. For that reason, none 
of our simulations reproduces the data in the upper right cor- 
ner of the figure. 

Figure [6] shows a comparison between our results and the 
observed correlation between the isotropic energy and the 
Lorentz factor of the outflow (Liang et al. 2010; Ghirlanda 
et al. 2012a). The observational correlation has significant 
scatter, possibly due to the difficulty of precisely measuring 
the Lorentz factor of the outflow. As an example of such dif- 
ficulty we show the measurements under the assumption of a 
wind environment with dark empty symbols and the measure- 
ments under the assumption of uniform environment with a 



Fig. 6. — Comparison between observed and simulated outflow Lorentz 
factor versus the burst energy. The black symbols with error bars show the 
measured Lorentz factor versus the isotropic equivalent energy, assuming that 
the external shock propagates into a wind environment. The colored lines 
and symbols show the result of our simulations, quantitatively reproducing 
the observed correlation. Colors and symbols are the same as in Figure |5| 
The Lorentz factors from the simulation are computed at the photosphere, 
since after the photons are released there is no internal pressure left to pro- 
duce further acceleration of the ejecta. The inset shows the Amati correlation 
highlighting in red the bursts for which the Lorentz factor could be measured. 
Finally, smaller gray symbols show the Lorentz factor derived from observa- 
tions under the assumption of a uniform external medium. 

solid gray symbol (from Ghirlanda et al. 2012). The two re- 
sults are indeed quite different. The Lorentz factor of the syn- 
thetic curves is computed as the average of the Lorentz factor 
of the material at the photospheric radius weighted with the 
local emissivity. Beyond the photosphere the flow is not ac- 
celerated since the radiation decouples and the outflow coasts 
until the external shock radius is reached, when the obser- 
vational Lorentz factor is measured. Our simulations repro- 
duce the correlation for a wind environment both in slope and 
normalization, but they have a smaller scatter. Simulated jets 
have, however, lower Lorentz factor than the measured val- 
ues for a uniform environment. Our results extend to lower 
Lorentz factors with respect to the measurements because it 
has not been possible to measure the Lorentz factor of weak 
GRBs, so far (see the inset in Figure [6]). The agreement of the 
simulations with the observations is entirely due to the jet-star 
interaction, since the simulations adopt a constant Foe = 400 
input except for one case, for which Too = 100 (the blue line 
and symbols). After the interaction with the stellar material 
has taken place, the Lorentz factor of the outflow at the pho- 
tosphere is strongly dependent on the polar angle, in such a 
way that the observed correlation is reproduced. 

Finally, we can use our simulated light curves to compute 
the radiative efficiency of the synthetic ORB prompt emission. 
We define the radiation efficiency as the ratio of the energy 
released as radiation at the photosphere over the total kinetic 
energy of the flow. We find (see Figure |7]) that the simula- 
tions predict a broad correlation between peak frequency and 
radiative efficiency and, transitively, between the total energy 
radiated and the efficiency. 

Measuring the local radiative efficiency from observations 
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Fig. 7. — Local radiative efficiency of the photospheric emission, defined 
as the ratio of the energy released in photons to the total energy of the out- 
flow. A tendency of bursts with lower peak frequency (and therefore with 
lower isotropic-equivalent energy) to have lower efficiency is observed in the 
simulations. Solid black symbols show local radiative efficiencies estimated 
by comparing the prompt and early afterglow emission of observed long- 
duration GRBs. 

is not an easy task (see, e.g., Zhang et al. 2007), since the 
measurement must be performed while the outflow is stiU rel- 
ativistic in order to avoid contributions from the ejecta moving 
at large angles with respect to the line of sight. For this reason, 
radio calorimetry cannot be used. The best proxy to the lo- 
cal radiative efficiency can be obtained, under some assump- 
tions, by comparing the brightness of the prompt emission 
with that of the early X-ray afterglow. Within the framework 
of an external shock dominated afterglow, the X-ray emission 
at early times is produced by electrons in the cooling regime 
and is fairly insensitive to the density of the ambient medium. 
Adopting the analytic approximations of Panaitescu & Kumar 
(2000) for a wind environment with electrons in slow cooling, 
we can derive the kinetic energy left in the fireball after the 
prompt emission as: 

where p is the slope of the non-thermal distribution of elec- 
trons, Ljj is the luminosity density at a frequency v above the 
cooling frequency, and eg are the shock equipartition pa- 
rameters (i.e., the fraction of the downstream internal energy 
used to accelerate relativistic electrons and to generate or am- 
plify the magnetic field, respectively), and Td is the time in 
days after the GRB trigger. Note that Eq. [3] depends very 
weakly on eg for a reasonable value of the electron acceler- 
ation slope p ^2 and only the electron equipartition param- 
eter plays a dominant role in the resulting calculation of the 
kinetic energy (Kumar & Piran 2000). 

In Figure'^ we show the observational efficiencies from a 
sample of GRBs with X-ray afterglows measurements with 
solid black symbols. We compute the £'x,iso released in the 



time interval (200,250) s, rest-frame using the best-fitting 
light-curves profiles obtained by Margutti et al. (2013) in the 
common rest-frame energy band 0.3-30 keV. Note that this 
estimate does not include the contribution from X-ray flares. 
Only GRBs with early afterglow re-pointing have been in- 
cluded and no extrapolation of the observed temporal behav- 
ior has been performed. The resulting efficiencies depend on 
the assumed value and therefore the agreement observed 
in Figure |7] between the data and observations should be 
seen as a consistency check rather than a proof. We adopted 
fairly standerd values in the calculations (e^ = 0.1, eg = 0.01, 
/7 = 2.3), but it has to be acknowledged that the quantitative 
agreement would be lost for a smaller or bigger value of (a 
qualitative good agreement is obtained for 0.05 < < 0.2). It 
is worth mentioning, though, that the general trend of higher 
efficiency for bursts with larger peak frequency is robust and 
does not depend on the choice of e^. Finally it should be 
stressed that at the considered time (200 - 250) s there might 
still be some contribution to the X-ray light curve form the 
prompt emission. We have also tried to use data from an ear- 
lier (100 < ^ < 150 s) and a later (500 < r < 750 s) interval, 
obtaining analogous results. 

5. DISCUSSION AND CONCLUSIONS 

We presented the results of a set of hydrodynamic numer- 
ical simulations of the evolution of long-duration GRB jets 
and their progenitor stars. Thanks to the extended domain of 
the simulations we were able to derive the location at which 
the spectral peak of the advected radiation is set and use that 
information to compute light curves from different progeni- 
tors and as seen by observers at different off-axis angles with 
respect to the stars' rotation axis. We find that within the 
assumptions made in this study, GRBs produced by photo- 
spheric emission reproduce some observational correlations 
such as the Amati correlation and the energy-Lorentz fac- 
tor correlation. The ability to account for such correlations, 
where other radiation mechanisms have failed, suggests that 
the bulk of the GRB prompt emission is released at the photo- 
sphere. Our results also suggest that the primary driver of the 
difference between very energetic events and weaker ones is 
the viewing angle, while the scatter in the Amati correlation 
is due to different jet/progenitor configurations. 

Our simulation and the radiation code applied in post- 
processing do not follow the full radiation transfer through 
the outflow and therefore the radiation properties were com- 
puted assuming that the spectrum of the entrained radiation 
is well described at all times by a Planck function in ther- 
mal equilibrium with the leptons. In a non-magnetized rel- 
ativistic jet photon production mechanisms are severely sup- 
pressed and the outflow becomes scattering-dominated well 
before the photospheric radius (Beloborodov 2013; Vurm et 
al. 2013). Under such conditions, especially in the presence 
of dissipation, the radiation spectrum would be better approxi- 
mated by a Wien function with a harder F{y)(xu^ low-energy 
tail. We note that, as long as the radiation and leptons are in 
thermal equilibrium, the peak frequencies of the photons in 
a Planck and Wien spectrum are equal to within a small fac- 
tor i^peak,Pianck/^peak,wien = 1-07. Our results are therefore only 
marginally affected by the choice of a Planck spectrum ver- 
sus a Wien one. We also note that a detailed discussion of 
the spectral shape of the photospheric emission is beyond the 
scope of this paper and we refer the reader to the extensive 
literature on the formation of non-thermal spectra when sub- 
photospheric dissipation and/or magnetic fields are present 
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(Giannios 2006; Pe'er et al. 2006; Giannios & Spruit 2007; 
Beloborodov 2010; Lazzati & Begelman 2010; Vurm et al. 
2011). 

Another result of our simulations is a strong correlation be- 
tween the average Lorentz factor of the ejecta that contribute 
to the light curve and the total radiated energy in the direc- 
tion to the observer. Such a correlation has been observed 
(Liang et al. 2010; Ghirlanda et al. 2012) and the results of 
our simulations are in good agreement with observations as 
long as a wind environment is assumed in the calculation of 
the observational Lorentz factor. This is another important 
piece of the puzzle, since the "standard" synchrotron internal 
shock model predicts the opposite correlation (Ghisellini et al. 
2010). Alternatively, the correlation might be related to the 
entrainment of baryons from a neutrino-driven wind (Lei et 
al. 2012). Finally, we found in our simulations that the radia- 
tive efficiency of brighter bursts is higher than that of weaker 
bursts. By comparing the early X-ray afterglow flux to the 
prompt emission energy of a sample of well-observed bursts, 
we were able to show that the correlation exists in the obser- 
vations (see also Bemardini et al. 2012). However, a detailed 
comparison requires the assumption of a value for the electron 
equipartition parameter e^, making the quantitative compari- 
son with the simulations parameter-dependent. We show that 
adopting the fiducial value = 0. 1, we can obtain good agree- 
ment between the simulation results and the observed values. 



As a final note it should be remembered that the correla- 
tions that we discuss here are not exhaustive of all correla- 
tions discovered in GRB catalogs. In particular, there is spec- 
ulation that the more physical correlation is between the peak 
frequency and the luminosity, rather than the energy release 
integrated over time (Yonetoku et al. 2004). In addition, an 
analogous correlation with luminosity exists within individ- 
ual bursts, where the instantaneous peak frequency correlates 
with the instantaneous luminosity (Lu et al. 2012). Since 
our simulations assume, for lack of a better model, an engine 
with constant luminosity, we cannot explore any correlation 
that depends on the variability injected by the central engine. 
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